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ABSTRACT 


We  present  in  this  study  a  three  step  algorithm  for  the  decomposition  of 
arbitrary,  three-dimensional,  planar  polygons  into  convex  polygons.  Through  a 
series  of  translations  and  rotations,  an  arbitrary  polygon  is  mapped  onto  the  x-y 
plane,  then  broken  into  a  set  of  convex  polygons,  and  finally  mapped  back  to  the 
polygon's  original  coordinate  system  for  filling  and  display  by  special  graphics 
hardware. 
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I.  INTRODUCTION 


Computer  hardware  advancements  applied  to  the  field  of  computer  graphics 
have  recentlv  made  possible  the  real-time  display  and  animation  of  three- 
dimensional  object  models.  In  most  currently  produced  graphics  systems,  the 
three-dimensional  (3D)  object  models  are  comprised  of  sets  of  3D  polygons.  For 
years  the  key  bottleneck  to  generating  real-time  images  of  3D  models  has  been  the 
limited  speed  with  which  the  3D  polygons  comprising  the  model  have  been  filled. 
This  bottleneck  has  been  alleviated  recently  with  the  advent  of  special  hardware 
designed  to  fill  polygons  at  rates  exceeding  44  million  pixels  per  second  [Ref.  1], 
For  real-time  animation  of  3D  object  models,  such  polygon  fill  rates  are  essential. 

The  majority  of  the  graphics  systems  that  provide  polygon  fill  hardware  only 
provide  hardware  for  filling  convex  polygons.  Concave  polygons  passed  to  such 
hardware  usually  are  filled  incorrectly,  if  at  all.  To  utilize  the  capabilities  of  such 
hardware  for  the  general  polygon  case,  that  is,  a  mix  of  concave  and  convex 
polygons,  it  is  necessary  to  have  some  means  of  decomposing  an  arbitrary  3D 
planar  polygon  into  a  set  of  3D  planar  convex  polygons  completely  covering  the 
area  defined  by  the  original  polygon.  The  following  study  attempts  to  provide 
such  an  algorithm.  We  then  discuss  the  implementation  and  performance 
characteristics  of  that  algorithm  on  the  Silicon  Graphics.  Inc.  IRIS  workstation. 


II.  BACKGROUND 


A.  DEFINITIONS 

A  polygon  is  defined  as  a  geometric  figure  that  is  described  by  an  array  of  x. 
v.  and  z  coordinates  for  its  vertices.  A  two  dimensional  polygon  lying  in  the  x-y 
plane  is  a  special  case  of  a  three  dimensional  polygon  in  which  th^  z  coordinates 
are  all  zero.  A  convex  polygon  is  one  that  has  no  points  of  negative  curvature 
along  its  boundary  [Ref.  2:  p.  217],  This  means  that  a  line  can  be  drawn  from 
any  point  inside  the  polygon  to  any  of  its  vertices  and  have  the  line  lie  entirely 
within  the  polygon  boundary  [Ref.  3:  p.  350).  A  concave  polygon  fails  the  above 
test.  i.e.  we  cannot  draw  a  line  from  every  point  inside  the  polygon  to  every 
vertex  and  have  the  line  remain  within  the  polygon  boundary.  However,  any 
point  inside  a  concave  polygon  can  be  shown  to  lie  inside  a  convex  polygon 
created  by  some  of  the  vertices  and  sides  of  the  original  polygon  [Ref.  3:  p.  330]. 
Figure  2.1  provides  examples  of  convex  and  concave  polygons. 

B.  DECOMPOSITION 

To  fill  a  concave  polygon  utilizing  hardware  that  only  correctly  fills  convex 
polygons,  we  must  first  break  down  the  polygon  into  simpler  convex  polygons  that 
together  cover  the  area  of  the  original.  The  process  of  decomposing  a  concave 
polygon  does  not  necessarily  yield  a  single  unique  set  of  convex  polygons,  but  the 


result  of  two  different  decompositions  of  the  same  polygon  should  be  equivalent  as 
shown  in  Figure  2.2. 

During  the  operation  of  decomposing  a  polygon,  there  are  two  main  checks 
that  must  be  conducted  to  ensure  complete  and  proper  decomposition  takes  place. 
The  first  check  is  to  ensure  that  no  area  of  a  created  convex  polygon  falls  outside 
the  boundary  of  the  original  concave  polygon.  The  second  check  is  to  ensure  that 
the  concave  polygon  is  completely  covered  by  the  total  areas  of  the  created  convex 
polygons.  Figure  2.3  provides  examples  of  the  two  checks. 

The  polygon  decomposition  algorithm  produced  for  this  study  is  comprised  of 
three  basic  steps.  The  first  maps  the  three  dimensional  coordinates  of  the  polygon 
into  a  two  dimensional  plane.  This  process  generates  a  polygon  that  is  more 
easily  examined  and  decomposed.  The  second  step  conducts  the  actual 
decomposition  an  performs  the  two  decomposition  checks.  The  third  step  maps 
the  coordinates  of  the  created  convex  polygons  back  to  the  the  original  orientation 
in  three-space.  At  this  point,  the  created  convex  polygons  are  then  ready  for 
passing  to  the  graphics  system's  polygon  fill  hardware. 

C.  RESEARCH  FACILITIES 

The  IRIS  (Integrated  Raster  Imaging  System)  2400  Graphics  Workstation, 
manufactured  by  Silicon  Graphics.  Inc.,  is  the  target  system  for  our  polygon 
decomposition  algorithm.  By  incorporating  custom  VLSI  chips  to  replace  slower 
graphics  software,  the  IRIS  svstem  can  accept  polygons  in  user-defined 


Figure  2.3a  -  Without  a  check  to  ensure  the  convex 
polygon  does  not  fall  outside  the  bounds  of  the 
original  polygon,  vertices  could  be  chosen  that 
yield  an  incorrect  decomposition. 


coordinates,  and  transform  them  to  screen  coordinates  for  display  on  a  color 
graphics  monitor  that  provides  1024  x  768  resolution.  This  use  of  special 
hardware  yields  processing  speeds  that  are  some  100  times  faster  than  similar 
operations  carried  out  in  software. 
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III.  STEP  I:  MAPPING  THE  3D  POLYGON  TO  THE  2D  X-Y  PLANE 

A.  THE  NORMAL  TO  A  POLYGON 

The  first  step  of  the  polygon  decomposition  algorithm,  is  comprised  of  several 
small  steps.  The  goal  of  these  steps  is  to  map  the  polygon  to  the  x-v  plane.  First, 
three  non-colinear  points  are  selected  from  the  array  of  x.  y.  and  z  coordinates  of 
the  polygon.  These  are  used  to  define  the  equation  of  the  plane  for  the  3D 
polygon,  and  the  normal  to  that  plane  (Figure  3.1). 

Ax  +  By  +  Cz  =  D 
N  =  Aj  +  Bj  +  Cg 

The  coefficients  A.  B.  and  C  are  calculated  from  the  coordinates  of  the  three 
selected  points  by  the  following  equations: 

A  =  (y0  -  yi)(*j  -  *i)  -  (ya  -  yi)(*o  -  *1): 

B  --  (z0  Z,)(x,  x,)  -  (z,  z,)(Xo  x,): 

C  =  (xo  x,)(yj  y,)  (xs  x,)(y0  y ,  >. 

The  coefficients  are  then  used  to  solve  for  D  of  the  equation  of  the  plane. 


B.  ANGLES  ABOUT  THE  AXES 

Once  the  equation  of  the  plane  is  determined,  each  point  in  the  polygon  input 
array  is  checked  to  ensure  that  the  polygon  is  in  fact  planar.  If  the  polygon 
passes  this  test,  then  the  angles  between  the  normal  and  the  axes  of  the 
coordinate  system  are  found.  From  these  angles,  the  necessary  rotations,  about 
the  x  and  y  axes,  to  bring  the  polygon  into  the  x-y  plane,  are  calculated.  The 
following  equations  [Ref.  4:  pp.  24-59]  are  used  to  calculate  the  angles  shown  in 
Figure  3.2a. 


cos  “  "  v/FTFTc5 
COS  *  =  y/A2  +  B*  +  C2 


COS  7  - 


_ C _ 

v/atTbtTc* 


For  the  polygon  to  lie  in  the  x-y  plane,  the  angle  between  the  normal  and  the 
z  axis  is  either  zero  or  180  degrees.  Necessarily  then  rotations  are  conducted  to 
bring  the  angles  between  the  normal  and  the  x  and  y  axes  to  90  degrees  each. 

First  a  projection  of  the  normal  onto  the  x-z  plane  is  made  as  shown  in  Figure 
3.2b.  The  dot  product  of  that  projected  normal  .V  and  the  z  axis  is  calculated. 
The  dot  product  is  divided  by  the  product  of  the  magnitudes  of  the  two  vectors 
yielding  the  cosine  of  the  angle  of  rotation  about  the  y  axis  (</>). 
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A  second  projection  of  the  normal  onto  the  y-z  plane  is  made  and  the  angle  of 
rotation  about  the  x  axis  (0)  is  calculated  in  the  same  manner  as  above  (Figure 
3.2c). 

C.  TRANSFORMATION  MATRIX 

To  actually  rotate  the  polygon  into  the  x-y  plane,  a  transformation  matrix 
must  be  created  that  is  used  in  the  multiplication  of  the  polygon's  coordinates. 
We  first  translate  the  polygon  to  the  origin  by  adding  the  negative  x.  v.  and  z 
coordinates  of  one  of  the  polygon's  vertices.  Assuming  1.  m.  and  n  are  the 
coordinates  selected,  the  matrix  becomes: 

10  0  0 
0  10  0 
0  0  10 

-I  m  -n  1 

A  rotation  matrix  for  the  angle  <t>  about  the  y  axis  is  then  created 


cos^  0  sin^  0 
0  10  0 
-sin^  0  cos^  0 
0  0  0  1 

and  concatenated  with  the  translation  matrix.  This  is  followed  by  a  multiplication 


with  a  rotation  matrix  for  the  angle  6  about  the  x  axis. 


10  0  0 
0  cos#  -sin#  0 
0  sin#  cos#  0 
0  0  0  1 

The  result  of  these  matrix  concatenations  is  the  following  transformation  matrix 
that  is  used  to  map  the  array  of  vertices  onto  the  x-v  plane. 


ros^ 

-  sin^sin# 

-sin^cos# 

0 

0 

cos# 

sin# 

0 

sin^ 

cos#sin# 

cos^cos# 

0 

lros^  nros# 

moos#  +  (lsin^-  ncos^) (  sin#) 

msin#+(lsin^  ncos#)(cos#| 

1 

Once  we  have  this  matrix  set  up.  call  it  M,oplant.  we  then  use  it  to  transform 
the  polygon's  coordinates  via  a  simple  row  vector  by  matrix  multiplication.  The 
result  of  this  multiplication  is  a  transformed  set  of  coordinates,  in  which  the  z 
coordinates  are  zero.  The  x  and  v  values  of  this  transformed  array  are  then  used 
in  tht  decomposition  step  of  our  algorithm.  We  note  at  this  point  that  for  the 
mapping  of  the  coordinates  back  to  the  original  vertices  (Step  III  of  the  Alg.j. 
that  the  index  into  the  transformed  array  is  also  the  index  for  the  untransformed 
point  in  the  original  array  of  vertices. 


IV.  STEP  II:  POLYGON  DECOMPOSITION 


The  decomposition  step  divides  the  polygon  into  multiple  convex  polygons. 
In  so  doing,  the  algorithm  performs  the  two  checks  to  ensure  that  no  area  outside 
the  polygon  is  shaded  and  that  the  created  convex  polygons  completely  cover  the 
area  of  the  concave  polygon.  This  step  requires  that  the  vertices  of  the  polygon 
be  ordered  in  a  clockwise  fashion,  with  respect  to  the  +z  directed  normal  to  the 
x-y  plane.  Thus  a  determination  of  the  ordering  must  be  made  and  if  necessary 
the  array  of  vertices  reversed. 

A  CLOCKWISE  ORDERING 

To  determine  if  the  vertices  of  the  polygon  are  ordered  clockwise,  the 
polygon's  area  is  calculated  using  the  formula  [Ref.  5:  p.  369]  below.  (Note  r,.y, 
are  the  polygon's  2D  coordinates.) 

yUx^y,)  t  (x,y,)  t  •••  t(x„  iy„)  Mx„y0)  (y0x,)  (y,x,)  -(y„-|Xj  (ynx<,)) 

The  above  equation  yields  a  negative  area  for  clockwise  and  a  positive  area  for 
counterclockwise  polygons.  If  the  polygon  is  ordered  counterclockwise  then  the 
array  of  x.  y.  and  z  coordinates  that  define  the  polygon  is  reversed  before  the 
selection  procedure  is  begun. 
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B.  VERTEX  EXAMINATION 

The  selection  of  vertices  for  a  convex  polygon  is  conducted  by  examining  each 
vertex  of  the  concave  polygon  with  respect  to  the  previously  selected  vertices  of 
the  developing  convex  polygon.  This  is  begun  by  initially  selecting  two  adjacent 
vertices,  then  in  clockwise  order,  examining  each  of  the  polygon's  remaining 
vertices  to  build  as  large  a  convex  polygon  as  will  fit  within  the  boundary  of  the 
concave  polygon. 

A  series  of  tests  are  conducted  in  the  selection  of  vertices  for  the  convex 
polygon.  As  vertices  are  selected,  the  line  segments  connecting  these  become  the 
edges  of  the  convex  polygon  and  define  its  boundary. 

1.  Test  One 

The  first  test  of  a  vertex  is  if  it  lies  to  the  left  or  right  of  each  line 
segment  of  the  developing  convex  polygon.  The  vertices  of  the  concave  polygon 
are  examined  in  clockwise  order.  Thus  if  a  vertex  lies  to  the  left  of  any  of  the  line 
segments  of  the  developing  polygon,  then  a  concave  angle  is  necessary  to  include 
the  vertex  in  the  polygon.  As  shown  in  Figure  4.1a.  vertex  2  lies  to  the  left  of  the 
line  segment  0-1.  A  concave  angle  (  p  in  Figure  4.1b)  is  necessary  to  add  vertex  2 
to  the  developing  polygon.  Therefore  vertex  2  is  not  selected  and  the  process 
moves  on  to  perform  the  same  test  on  vertex  3. 

2.  Test  Two 

A  vertex  that  passes  the  first  test  may  still  require  a  concave  angle  to  be 
included  in  the  developing  convex  polygon.  In  Figure  4.2.  the  vertices  0.  1.  2.  3. 


Figure  4.1  -  (a)  Vertex  2  lies  to  the  left  of  the 
line  segment  0-1,  (b)  Concave  angle  ( p )  would 

be  required  to  connect  Vertex  2  to  the  devel¬ 
oping  polygon,  (c)  Vertex  2  is  skipped  and 
the  examination  is  continued  with  vertex  3. 
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Figure  4.2  -  (a)  Vertex  5  lies  to  the  right  of  each 
line  segment  of  the  developing  polygon,  (b)  Concave 
angle  ( p )  is  required  to  connect  Vertex  5  to  Vertex 
0,  (c)  Concave  angle  ( p )  can  be  identified  since  0 

lies  to  the  left  of  line  segment  4-5,  (d)  Vertex  4 

is  skipped  and  the  test  repeated  with  line  segment 
3-5,  (e)  Vertex  3  is  skipped  and  the  test  repeated 

with  line  segment  2-5,  (f)  Result  of  the  tests  are 

that  vertices  0,  1,2  and  5  are  selected. 


and  4  have  been  previously  selected  for  a  convex  polygon  with  vertex  5  under 
examination.  Vertex  5  passes  the  first  test  since  it  lies  to  the  right,  or  inside  of 
each  of  the  line  segments.  However,  a  concave  angle  (f>  in  Figure  4.2b)  is  again 
required  to  connect  vertex  5  with  vertex  0  of  the  developing  polygon.  A  concave 
angle  of  this  type  is  identified,  as  shown  in  Figure  4.2c.  if  the  initial  vertex  of  the 
developing  convex  polygon  lies  to  the  left  of  the  line  segment  created  by  the 
vertex  under  examination  and  the  last  selected  vertex  of  the  developing  polygon. 
In  Figure  4.2c.  the  initial  vertex.  0.  lies  to  the  left  of  line  segment  4-o.  To  remedv 
this,  the  last  selected  vertex  of  the  developing  convex  polygon,  vertex  4.  is 
removed  and  the  test  is  again  conducted  with  the  line  segment  formed  by  the 
vertex  under  examination  and  the  last  selected  vertex  (Figure  4. 2d  4.2e).  This 
process  is  repeated,  removing  vertices  from  the  developing  polygon,  until  the 
initial  vertex  no  longer  lies  to  the  left  of  a  new  line  segment. 

3.  Test  Three 

The  third  test  conducted  is  to  ensure  that  no  line  segment  created  in  the 
convex  polygon  intersects  any  of  the  edges  of  the  concave  polygon.  Such  an 
occurrence  results  in  area  outside  the  concave  polygon  being  included  in  the 
convex  polygon  and  an  inaccurate  polygon  fill  taking  place.  In  Figure  4.3.  the 
previously  selected  vertices  are  0.  1.  2.  and  3.  Vertex  7  is  under  examination  and 
passes  the  two  earlier  tests.  However,  the  -election  of  vertex  7  results  in  the 
convex  polygon  0.  1.  2.  3.  7  which  contains  area  outside  the  concave  polygon.  A 
situation  such  as  this  is  identified  since  line  segment  3-7  intersects  line  segment 
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Figure  4.3  -  Tests  must  be  conducted  for  the 
intersection  of  developed  line  segments  with 
the  edges  of  the  original  concave  polygon. 
Otherwise  Vertex  7  would  be  an  acceptable 
vertex  for  inclusion  in  the  developing 
polygon  0 ,  1 ,  2 ,  3 . . . 


13-0  of  the  concave  polygon.  When  this  occurs,  the  vertex  under  examination  is 
removed  from  the  convex  polygon  and  the  next  vertex  of  the  concave  polygon  is 
examined. 

(’.  SELECTION  COMPLETION 

The  selection  process  for  the  developing  convex  polygon  concludes  when  all 
vertices  of  the  concave  polygon  have  been  examined.  The  convex  polygon  is  then 
ready  for  the  third  step  of  our  algorithm,  mapping  to  the  original  polygon's  3D 
coordinates.  Figure  4.4  is  an  example  of  a  concave  polygon  that  was  decomposed 
to  the  three  shaded  convex  polygons  using  the  above  three  tests.  We  have  a 
"hole"  left  in  the  concave  polygon  even  though  all  vertices  and  edges  of  the 
concave  polygon  have  been  utilized.  Holes  such  as  this  are  prevented  by  taking 
the  initial  and  last  vertices  of  a  developed  convex  polygon  as  the  starting  vertices 
of  a  convex  polygon  to  be  developed  In  Figure  4.4.  vertices  0  and  2  of  polygon  0. 
1.  2  are  taken  to  develop  convex  polygon  0.  2.  4  to  fill  in  the  "hole".  Once  this  is 
accomplished,  step  two  of  the  algorithm  is  repeated  by  selecting  a  vertex  of  the 
concave  polygon  that  has  not  yet  been  utilized  in  a  convex  polygon.  The 
decomposition  process  is  complete  when  there  are  no  vertices  or  edges  of  the 
concave  polygon  that  have  not  been  utilized  in  one  or  more  convex  polygons. 
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V.  STEP  III:  MAPPING  TO  THE  ORIGINAL  COORDINATE  SYSTEM 

The  second  and  third  steps  of  the  decomposition  process  are  interleaved,  i.e.. 
each  created  convex  polygon  is  immediately  mapped  back  to  the  original 
coordinate  system.  This  is  accomplished  by  maintaining  the  original  array  of 
vertices  that  the  user  provided  for  the  three-dimensional  polygon.  As  a  convex 
polygon  is  produced,  it  is  actually  represented  by  the  indices  of  its  vertices.  Actual 
vertex  coordinates  are  extracted  from  the  vertex  array,  as  required,  using  the 
index.  In  this  manner,  the  transformed  coordinates  can  be  utilized  during 
decomposition  and  the  index  can  be  used  to  obtain  the  original  coordinates  for 
final  display.  By  this  method,  there  is  no  consequential  need  to  translate  and 
rotate  the  vertices  back  to  the  original  coordinate  system. 


VI.  IRIS  IMPLEMENTATION 

Implementation  of  the  decomposition  algorithm  on  the  IRIS  system  is  via 
approximately  1600  lines  of  C.  code  (Appendix).  Several  data  structures  and 
various  IRIS  system  hardware  features  are  utilized  as  described  below. 

A.  MATRIX  IMPLEMENTATION 

In  the  development  of  the  transformation  matrix.  IRIS  system  calls  for 
translation  and  rotation  are  used  to  alter  the  top  of  the  system  matrix  stack. 
These  routines  perform  the  matrix  concatenations  needed  to  translate  and  rotate 
the  polygon  to  the  x-y  plane.  To  convert  the  polygon  vertices  to  2D  coordinates, 
four  vertices  at  a  time  are  loaded  into  a  matrix.  The  x,  y,  and  z  coordinates  go 
into  the  first  three  columns  with  l's  filling  the  fourth  column  to  provide  the 
necessary  homogeneous  coordinates.  This  is  placed  on  the  system's  matrix  stack 
and  multiplied  with  the  transformation  matrix  using  the  IRlS's  Geometry 
Engines.  The  transformed  coordinates  are  removed  from  the  matrix  stack  and  the 
process  is  repeated  until  all  vertices  have  been  transformed 

B  DOI'BLY  LINKED  LIST 

A  data  structure  for  each  vertex  of  the  polygon  is  maintained  in  our  system. 
This  structure  contains  the  vertex’s  index  into  the  vertex  array.  As  actual 
coordinates  are  needed,  the  vertex  s  index  is  used  to  extract  the  x.  y.  and  i. 
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coordinate*  from  this  array.  A  "used"  flag  indicates  the  vertex's  use  in  a 
developed  convex  polygon,  with  a  slope  and  y-intercept  stored  for  the  line 
segment  the  vertex  forms  with  its  predecessor  in  the  array.  A  forward  pointer 
connects  the  structure  with  the  next  vertex  in  clockwise  order,  while  a  back 
pointer  connects  the  structure  with  its  preceding  vertex.  This  produces  a  doubly 
linked  circular  list  of  vertex  structures  that  permits  forward  clockwise  traversal  of 
the  polygon  and  rearward  counterclockwise  traversal  (Figure  6.1). 

(  .  VERTICES  STACK 

To  manage  the  examination  and  selection  of  vertices  for  a  developing  convex 
polygon,  a  stack  is  utilized.  The  size  of  the  vertices  stack  is  determined  at  run 
time  when  contiguous  memory  is  dynamically  allocated,  based  on  the  number  of 
vertices  in  the  concave  polygon.  As  a  vertex  is  examined,  its  index  is  placed  on 
top  of  the  stack  along  with  the  slope  and  y-int.ercept  of  the  line  segment  it  creates 
with  the  vertex  below  it.  If  the  vertex  passes  all  selection  tests,  it  remains  on  the 
stack  and  a  new  vertex  is  pushed  on  top.  On  the  other  hand,  when  a  vertex  fails 
a  selection  test,  it  is  popped  off  the  stack  and  the  next  vertex  to  be  examined  is 
placed  on  top.  Once  the  selection  of  vertices  for  a  convex  polygon  is  completed, 
the  vertices  are  passed  to  the  IRIS  system  hardware  and  the  stack  is  cleared  to 
continue  with  the  decomposition  process. 
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VII.  LIMITATIONS  ANDC ONC L U S ION S 

This  study  has  presented  an  algorithm  for  the  decomposition  of  an  arbitrary, 
planar  polygon  into  a  set  of  convex  polygons,  thus  eliminating  the  need  to 
manually  perform  such  decomposition.  With  the  algorithm's  dynamic  allocation 
of  memory,  the  size  polygon  that  can  be  filled  is  limited  only  by  the  amount  of 
memory  available  for  use.  The  algorithm  indirectly  provides  an  extended 
capability  over  the  IRIS  system  s  limit  of  384  vertices  in  one  polygon  (Ref.  6:  p. 
2). 

However  due  to  the  time  required  to  test  each  vertex  for  concavity,  the 
algorithm  is  not  suggested  for  use  where  real-time  response  is  important.  To 
demonstrate  the  overhead  involved  in  testing  for  concavity,  a  convex  polygon  was 
filled  and  displayed  500  times  by  the  algorithm.  The  same  polygon  was  then 
filled  utilizing  only  the  IRIS's  polygon  fill  hardware.  The  results  are  provided  in 
Figure  7.1  for  various  polygon  sizes. 

Due  to  the  three  primary  steps  necessary  for  decomposition  and  the  tests 
conducted  on  each  vertex,  the  speed  of  the  algorithm  depends  entirely  upon  the 
number  of  vertices  and  concave  angles  of  the  polygon.  As  such,  the  success  of  the 
algorithm's  use  in  real-time  display  and  animation  is  based  upon  a  balance 
between  the  complexity  of  the  concave  polygons  used  and  the  speed  requirements. 
We  suggest  that  the  created  concave  polygons  be  computed  once  and  saved. 
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For  IRIS  usage,  this  is  perhaps  best  accomplished  through  the  use  of  the  "object 
mechanism  of  the  IRIS  graphics  library. 
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APPENDIX  ALGORITHM  SOURCE  CODE 


******»*#*:*:*:**.**;***#*************.*±*:*****X***********4.******** 

globals.h 

»*•.*•*•****»****«*.*.  «***»»***-vx**«***x*X****.*.«**.*»xxs.**4».**a:***.» 

*  include  ^'gl.h 

=  include  -  device. h  > 

=  include  >'stdio.h  ■ 

=  include  ■  rnath.h  ■ 

GLOBAL  CONSTANTS 

«««*4»*****»***«V***VX»-*4XX**X*X*:***X*X*XX*«.»*.*X**«A***»*«*4**» 

*  define  NT  LL  0 

r  define  BOTTOM  0 

=■  define  new(x)  ix  ‘)  malloc| sizeof( x) ) 

*  *******  *  *  *  •  •  •  ••*  *  ******  *********************************  *4**  • 

GLOBAL  TYPES 


typedef  int  boolean: 

’  Define  a  structure  to  contain  data  about  a  vertex  of  the  polygon 
These  structures  will  f>e  used  to  create  a  linked  list.  ‘ 

t  \  pedef  struct  point  { 
struct  point  ‘back; 
int  index; 
float  slope: 

Coord  intercept: 
int  direction: 
boolean  failed: 
boolean  segment  used: 
struct  point  ‘forward: 

I  . 

(• 

typedef  struct  point  point: 
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*  Define  a  structure  for  storing  data  about  each  vertex.  This  will 
be  used  to  create  a  stack  for  later  use.  ' 

typedef  struct  stack  { 
int  index; 
float  slope; 

Coord  intercept; 
int  direction. 

}: 

typedef  struct  stack  Stack; 

.  *  *  *  *  **********************************  ********  *  *  4  *4  *  *  *  X  w  *  *  *  *  4  4 

GLOBAL  DATA  STRUCTURES  AND  V  ARIABLES 

*44  •***«  4*4*4  *  *4*4  *4*  4  *******444  4  **************  *  *  *  *  4  *********** 

Stack  'stack: 

point  ’start  ptr.  "current  ptr: 
int  TOFU 

boolean  clockwise  error  flag; 
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ROUTINE:  convert. c 

««*m*******’t**4**:******:«:*».*'«:*:*.«.:a:«''xxx.x»xx  *  *«•.**»  *  »  ft  a  *  *  *  »  as 

This  routine  takes  th  initial  array  of  vertices,  passed  by 
the  user,  and  makes  the  calls  to  other  routines  to  perform  the 
necessary  translation  and  rotations  to  transform  the  array  into 
the  x-y-0  plane  This  new  array  is  then  used  bv  the  decomoser  to 
breakdown  the  polygon  into  convex  polygons. 


cpolff total  pts.  pts) 

int  total  pts: 

Coord  pts  3  : 

f 

'1 

int  pt  1.  pt2.  pt3.  cnt: 

float  A.  B.  C.  D.  anglex.  angley .  anglez.  'new  pts: 
Matrix  transform  1: 

Matrix  points: 

s'ack  (Stack  * )  calloc(total  pis  -  1.  sizeof(Stack)): 


dynamically  allocate  memory  to  create  an  array  to  hold  the 
vertices  of  the  newly  transformed  polygon 


new  pts  (Hoat  *)calloc  (3‘total  pts.  sizeof( float ) ): 


get  3  non-colinear  vertices  from  the  polygon  to  be  used  for 
calculating  the  normal  to  the  polygon. 


extract  3  pts(tota)  pts,  pts.  Apt  1 .  Apt2.  Apt3): 


calculate  the  coefficients  for  the  equations  to  the  plane 
of  the  polygon  and  the  normal  to  that  plane. 


determine  coeffs(total  pts.  pts.  pt  1 .  pt2.  pt 3 .  A  A.  A  ll.  AC.  AD) 


check  to  ensure  all  vertices  of  the  polygon  lie  in  the  same 
plane.  If  so  then  calculate  the  necessary  angles  of  rotation 

iff  planar  polygonf  total  pts.  pts.  A.  B.  C.  L> ) )  { 

calculate  the  rotation  about  the  y  axis  necessary  to 
bring  the  normal  into  the  y-z  plane. 

****»*»*-*-*******  +  ***i**im***.****.».»***»4a*  ******.**.■».» 

angles  from  normalfA.  B.  C.  ’>’.  (Scanglex.  Bangley): 

create  the  transformation  matrix  necessary  to  multiply 
with  the  array  of  vertices  to  generate  the  rotation. 

*•*  ***  +  »*4  +  *t»  +  ****■»*  »■*.**■••  **.4  4*xwa.M*4r-M.*.M****'*x*: 

build  transform  matrix! transform!,  angley,  y  .  pts.  ptl) 


perform  the  matrix  multiplication  with  the  array  of 
vertices  and  store  the  new  coordinates  in  the  new  array 
new  pts. 

map  ptsf total  pts.  pts.  transform  1.  points,  new  pts); 

now  extract  three  non-colinear  vertices  from  the  new 
array  of  vertices,  new  pts. 

•  ••t**«*««*»«A«4***»*X*  +  *****4--*4-i*«*2»:***».**4r**v».** 

extract  II  ptsftotal  pts.  new  pts. &pt  I.  &pt2.  &pt3): 

since  the  polygon  has  been  rotated  about  the  \  axis  we 
calculate  the  coefficients  for  the  equations  to  the 
polygon  one  more  time. 


calculate  the  rotation  about  the  x  axis 


•  s«ma*m****«.  *»»**»»»:*#»*%***i*«*********»***»»*»-*  ***.*»%.»*»»»«»».» 

angles  from  norma!(A.  B.  C.  'x‘.  Aanglex.  Aangley  I : 

create  the  transformation  matrix 


build  transform  matrix(transforml.  anglex. ‘x  .  new  pts.  ptl): 


perform  the  matrix  multiplication 


map  ptsltota!  pts.  new  pts.  transform!,  points,  new  pts): 


call  the  routine  which  will  perform  the  decomposition 
upon  the  newly  transformed  polygon  which  now  lies  in  the 
\-v-0  plane. (thus  to  he  handled  as  a  21)  polygon) 


split  (total  pts.  new  pts.  pts): 


once  t  he  dynamically  allocated  memory  for  the  array 
new-pts  has  served  its  purpose,  the  memory  is  freed. 


cfree) new  pts.  ({’total  pis.  sizeof( float ) ); 
rfreejstack.  total  pts.  sizeof(Stack )  |: 


if  riot  planar  t  hen  t  lie  decomposition  is  not  attempted. 


ret  urn. 


ROL  TINE:  extract. r 


r^***)*#************:*****^****  *  «  *  *  *  •  »  «  ••  *  *  •  ******  *  *»  *  •  *  *  ▼  »* 

This  routine  searches  the  array  of  vertices  to  find  three  points  which 

do  not  lit  the  same  line  and  returns  the  indices  of  these  point1- 
*»**♦**»*  ,.*:****.*  ***#**»**»**•:•»••*•»•****»<  «*••*•*****•• 


extract  3  ptsftotal  pts.  pts.  ptl.  pt2.  pt3i 

int  total  pts; 

<  ’oord  pts  ,'J  ; 

int  'ptl.  '  pt  2.  '  pt.'f: 

) 

t 

int  count .  finished; 

'ptl  (I;  ‘pi  2  (t;  '  pt.'f  0;  count  0:  finished  t); 

'  The  first  vertex  in  the  array  (index  ())  is  selected  for  the  first  point, 
whilelcount  total  pts  K  A  ffinished)  | 

'  VS  e  then  search  the-  rest  of  the  array  of  vertices  to  locate  two  more 
count  ■  I . 


iff  >2  0)  { 

’  If  the  second  point  has  not  been  selected  yet  then  we  check  that 
\.  y.  and  z  coordinates,  of  the  vertex  currently  under  exam,  are 
different. 


f(pts  ' 

pt  1  (1  ' 

L>ts  count  0 

pts  ' 

ptl  1  ' 

pts  count  1 

pts  ‘ 

ptl  2  ! 

pts  count  2  ) 

'  If 

so  then  the  vertex  becomes  the 

'  pt  2 

count; 

t 

f 

else  iff* pt 3  0)  { 

'  If  the  third  point  has  not  yet  been  selected,  then  we  check  that 
the  coordinates  of  the  vertex  under  examination  are  different 
from  both  the  first  and  second  point.  ' 

iff  (pts  *  pt  1  0  '  pts  count  0 
pt s  'ptl  I  '  pts  count  1 


39 


ptsi*ptl  2.  !=  pts  count  2  )  && 

(pts  *pt2i  0  !  =  pts  count  0 
pts  *pt2  1  !=  pts  count  1 
pts  *pt2  2  !=■  pts count  2  )) 

'*  If  the  coordinates  are  different  then  the  \ertex  is  selected 
as  the  third  point.  ’ 

*pt3—  count:  } 

else 

finished  I: 


ROl  TINE:  is  planar. c 


r*«4C**4«4******»***W***«4****4*4*44»**44*»*4*4«*44**4»«*J 


This  routine  takes  each  vertex  of  the  polygon  and  uses  it  to  solve  the 
equation  for  the  plane  as  calculated  earlier.  If  every  vertex  lies 
in  the  plane,  then  a  value  of  1  is  returned. 


>  «  a  «  4  •  *••«*****«**«******••***•«••*****  •  ***** 


planar  polygon! total  pts.  pts.  A.  B.  C.  I)) 


irit  total  pts: 
Coord  pts  : 
float  A.  B.  C.  I); 


int  count,  planar: 


count  0:  planar  1 : 


whilefcount  -  total  pts  A- A-  planar) 


The  equation  for  the  plane  is  setup  to  equal  zero  if  the  coordinates 
of  the  vertex  lies  in  the  plane.  A  range  of  .02  on  either  size  of 
zero  is  allowed  to  handle  round  off  error.  This  test  is  done  for 
each  vertex  in  the  polygon.  If  any  vertex  fails  the  test  then  the 
function  returns  (0)  for  non-planar.  ' 


if((A‘(pts  count  ())•  B*(pts  count  I  }•  C*(pts  count  2  )-  D)  .02 

(  A  *  ( pts  count  0  )  -  B'(pts  count  i  )*  C*(pts  count  2  )-  D)  -0i 
planar  0: 


count  •  I : 


return(planar) 


This  routine  uses  the  three  vertices  selected  from  the  polvgon  to 
calculate  the  coefficients  that  will  be  used  in  both  the  equation 
of  tne  plane  the  polvgon  lies  in  and  the  normal  to  that  plane. 


determine  coeffs(iotal  pts.  pis.  pt  1 .  pt'2.  pt3.  A.  B.  C\  D \ 

i  11 1  total  pts: 

Coord  pts  3  : 
int  pt  1.  pt2.  pt 3: 
tic,  it  *  A.  *  B.  *t  .  ■  I): 


tioat 

\U.  \1. 

\2.  >  0.  >  1 .  \  2.  zl).  z  1 .  z2:: 

\it 

pts  pt  1 

0  . 

>0 

pts  pt  1 

1  : 

zO 

pts  pt  1 

•>  ■ 

\1 

pts  pt2 

0  : 

>  1 

pts  pt 2 

1  : 

/I 

pts  pt 2 

•>  . 

\2 

pts  pt3 

t)  : 

>2 

pts  pt 3 

1  ; 

z  2 

pts  pt3 

2 

*A 

l>0->  1  1 

‘(z2-zl )  -  (>2->  1  )*(zO-zl ' 

‘If 

1  ZU-Z  1  ) 

1  x2-x  1 )  -  (z2-zl  I’lxO-xl 

‘C 

’  l > 2->  1  >  -  |\2-\U*f>0-\  1 

‘I) 

:  *  A  1  *  pts  pt  1  0  -  ("Bl’ptsptl 

rrt  urn; 

*(■  ‘pts  pt  l  2 


This  routine  utilizes  the  coefficients,  determined  from  the  normal  to 
the  polygon,  to  determine  the  angles  of  rotation  necessary  about  the  x 
and  y  axes.  These  rotations  are  necessary  to  bring  the  polygon  into  the 
x-y  plane.  This  routine  is  called  the  first  time  to  calculate  the 
rotation  about  the  v  axis.  The  second  time  it  is  called  is  to  determine 
the  x  axis  rotation.  The  parameter  "axis"  is  used  to  pass  which  angle 
is  to  be  calculated. 


angles  from  normal!  A,  B.  C.  axis,  anglex.  angley) 

float  A.  B. 
char  axis: 

float  ‘anglex.  'angley: 

float  pi.  degx.  degy.  degz: 

pi  3.1415926536: 

degx  acos(A  sqrtf  A*  A  •  B’ B  *  C*C))*360' (2’ pi): 

degy  acosfB  sqrtfA’A- B*B~C*0))*360  ( 2  *  pi ) : 

degz  acosfC  sqrtfA’  A  ~B*B  -  C"C))*360  f‘2'pi): 

iff  A  0  A  A  B  0)  { 

'  simply  translate  the  polygon  to  the  origin  * 

‘ anglex-  0.0: 

'angley  0.0:  } 

else  iff  A  0MC-  0)| 

rotate  polygon  90  degrees  about  x  axis  * 

'anglex  90.0: 

'angley  0.0:  } 

else  iff  B  0  A  A  O  0)  f 

'  rotate  polygon  90  degrees  about  v  axis  * 

'anglex  0.0: 

‘angle)  90.0:  } 

else  { 

must  calculate  amount  of  rotation  about  the  x  A  y  axes  ' 
‘anglex  acosfC  sqrt  (B*  B  *  C*C))  *360  ( 2  *  pi ) : 

‘angle>  acosfC  sqrtfA' A  •  C’*C'))’360  (2‘ pi); 


‘angie\  -  -( ’angle)  ): 
’  anglex  _  0.0:  } 
else  •{ 

if(deg>  90.0)  ■' 

' anglex --  -j  ’anglex  1: 
’angle)  -  0.0:  f 

v 


convert  the  angles  into  tenths  of  degrees  as  needed  bv  the  IR1.S 
‘anglex  10’ i  ‘anglex  :: 

*  angle)  10‘;  ‘angle)  : 


ROUTINE:  transmatx.c 

**********  4.  **1tn**.XX+**#X*  +  *  +  **  +  *4*1***  ***  +  tMtt  a 

This  routine  takes  a  matrix  structure  and  initializes  it  to  an  identity 
matrix.  Then  the  necessary  translation  values  and  rotation  values  are 
added  to  the  matrix  to  create  the  resultant  transformation  matrix. 

*«*••***««***  »•*•*»*:•***«««***»***••**  X**«»*«*  4J*  *«*«**«*• 

build  transform  rnatrix(transforml.  angle,  axis.  pts.  pt  1 ) 

Matrix  transform  1: 
float  angle: 
char  axis: 

Coord  pts  3  : 
int  ptl: 

i nt  i.  j: 

create  the  identity  matrix  ' 
fori  i  0:  i ■  4 :  i  -  I )  { 
for(j  0:  j-  1:  j  •  I  I  [ 

iff'-  j) 

transform  1  i  j  1: 
else 

transform  1  i  j  I): 

I' 

} 

pushmatrix( ): 
loadtnatrix(transforml ): 

add  in  the  rotation  values  * 
rotate((int  (angle,  axis): 

add  in  the  translation  values  ’ 
translate(-pts  pt  1  0  .  -pts  ptl  1. -pts  ptl  2); 

getmatrixf  transform  1 ): 
popmatrix( ): 

ret  urn: 
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HOI  TIN L:  mapover.c 


This  routine  places  up  to  4  vertices  into  a  matrix  at  one  1 1  me  then 
performs  a  matrix  muit ipitcalion  with  the  transformation  matrix  to 
produce  x.  v.  and  i  coordinates  for  the  polxgon  that  is  translated 
and  rotated  some  amount  about  the  x  and  v  axes  These  new  coordinate- 
are  restored  in  the  arrav  new  pts  for  further  manipulation. 


map  p;  total  pts.  pts.  transform!,  points.  new  pts| 

mi  total  pt - : 

(  oord  pts  d 

Matrix  transform!,  points. 

float  new  pts  . 


.tit  t  outit .  coutit'i.  flit : 
count  0:  count.’  0: 
w  hliei  count  tot  al  pts  •  \ 
i  nt  It; 

w  r.nel  (tit  1 

‘  mate  the  coordinates  of  a  vertex  into  the  matrix 
points  i  nt  (i  pts  count  0 

points  cm  1  pis  count  I  . 

point s  c nt  pt s  (  ount  '2 

points-  ent  .'f  1  0: 

if  I .  count  -  1  i  total  pts  • 

(Til  -  1  . 

eis< 

PreaK : 


save  the  state  of  the*  matrix  statk  and  perform  the  multiplication 
pushmatnxi  l. 

.oadrtiat  rixi  t  ransform  1  I , 
mult  mat  r:x  I  points  !: 

Set  mat  r t x t  points 


poprnatrixj ): 


rnt  0: 

while! nil  4)  { 

‘  move  the  newly  calculated  coordinates  back  into  the  vertices  array 


new 

pt  s  count  2 

0 

points 

cnt  0 

pt s  count  2 

1 

points  cnt  1 

new 

pt s  count 2 

u 

points 

cnt  2 

if(  ( count  2  1|  count ) 

cnt  •  1 : 

else 
break: 


ret  urn: 
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T 


t»l'l 


»r»i  ■» ■■ 


T 


HOI  TINK  split  c 


This  routine  tak'  -  the  arras  for  the  t  ransfortned  polygon  and  cai  l> 
the  t.eeesvars  fi.  .ions  and  routines  to  decompose  the  polygon 
and  pass  the  vrr  res  of  the  convex  polygons  to  the  I  R  1  >  system 
tor  tilling 


split  l  iota,  pt pt  s.  oat  pt  s  ! 


int  total  i>ts: 
(,'oord  pts  .'i  ; 

<  oord  old  pt  s  .'1 


Object  iiecoiiii>ose  po.vgom 
Object  whole  thing: 
boolean  c  lot  kw  ;se. 

Point  "pi.  '  pi!. 


Creat*'  the  doubly  ImketJ  circular  list,  "main  chain",  for  the 
concave  polygon  represented  by  the  pts  array.  \iso  initialize 
hag'  and  ptrs  that  will  t>e  used  in  decomposing  the  polygon. 


ciocswise  rlixkwise  ordering!  iota.  pts.  pt- 
create  rnait:  c  r  airu  ciockw  ise.  total  pt'.  pts  . 


Recompose  the  main  polygon  and  build  the  set  of  simple  convex 
; *c > I  v  goris 


de<  orr.fHise  poiv  gonl  total  pts.  pts.  i>ld  pts  : 
ret  urt. 
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\  v  -VP  -N 


*»****»***»r*****»:*:*********:r***.r4t  »*»*.****♦»  ******************* 

ROUTINE:  main  chain. c 

•.*****¥******#***  ****«C***»'*#*********^**»'******»**#**»*******» 

Create  the  doubly  linked  circular  list,  "main  chain",  for  the 
concave  polygon  represented  by  the  pts  array.  Also  initialize 
flags  and  ptrs  that  will  be  used  in  decompi  ,ing  the  polygon. 

«  *  *  *  »  *  •  *  »  4  «.*  *  ***»  <*»»*#»*•**♦»■***  **  »  *  «  »  •  s  a 


create  main  chain(clockwise.  total  pts.  pts) 

boolean  clockwise, 
ini  total  pts: 

( 'oord  pts  3  : 

I 


int  count  1 : 
point  ‘pi,  ‘p‘J: 

"  allocate  memory  for  a  structure  to  contain  data  on  the  first 
vertes  of  the  polygon  * 

pi  new  (point): 

'  set  pointers  to  this  structure  and  initialize  all  values  * 

start  ptr  pi: 

current  ptr  pi: 

pi-  index  0: 

pi-  direction  0: 

pi-  -failed  FALSE: 

pi-  segment  used  -  FALSE: 

'  if  polygon  vertices  are  ordered  clockwise  then  we  index  the  vertices 
array  beginning  at  0  and  incrementing  from  there.  Other  we  must  begin 
at  the  other  end  of  the  vertices  array  to  extract  our  coordinates.  In 
this  way  a  forward  move  along  our  doubly  linked  list  will  equal  a 
clockwise  move  around  the  polygon. 

iffclock  wise) 

count  I  1 : 
else 

( ount  1  total  pts-  1 : 

*  for  each  vertex  of  the  polygon  we  build  a  point  structure,  calculate 
the  slope,  intercept,  initilize  Hags  and  connect  the  structure 
via  pointers  to  its  preceeding  and  succeeding  vertices  in  the 
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polygon.  The  result  is  a  doublx  linked  list  of  struct ure'- 


do- 

p'J  new  (point); 
pi-  forward  -  p'2; 
p2-  hack  pi. 
pl  p'J: 

pi-  index  count  1: 

determine  slopeiptspl-  hack-  index  0  .  pts  pl-  ■  back-  index  1 
pis  pl-  index  0  .  pts  pl-  -index  1  .  Ajpl-  slope;. 

Alpl-  -intercept ).  A.' ( pl-  directionii: 
pl-  -failed  F  \ l.SK 
pl-  segment  used  FALSK: 
iflruK’kw  ise  ; 

count  1  -  1 

e.s< 

count  I  -  1 . 

px  niiei count  1  totai  pt»  A  A.  count  1  0:; 

‘  connect  :nc  nrst  xertex  structure  to  the  last  xertex  of  the  polxgon 
and  calculate  tne  slope  arid  intercept  of  the  iine  segment  the  two 
create . 

pl-  forward  start  ptr: 
star;  ptr-  hack  -  pl. 
pl  'tart  ptr. 

determine  slopeiptspl-  hack-  index  0  .  pts  pl-  hack-  .ndex  1  . 
pts  pl-  -index  (t  .  pts  pl-  index  1  .  k  I  pl-  -slope'-. 

A,pl-  intercept  i.  Ai  pl-  direction!!: 

ri  t  urn. 


5i 


$ 


KOI  TINE: 


clockwis.r 


f*-H***^-***4J 


This  routine  calculates  the  area  of  the  arbitrary  planar 
polygon  once  it  has  been  transformed  to  the  2D  x-y  plane. 

A  negative  area  indicates  a  clockwise  ordering  of  the  vertices 
while  a  positive  area  indicatesa  a  negative  orderig  of  vertices 
A  "0"  is  returned  for  counterclockwise  ordering  while  a  "1"  is 
returned  for  clockwise  ordering. 


clockwise  ordering) total  pts.  pts) 

int  total  pts: 

<  'oord  pts  .'1  : 


float  area: 
int  i.  result : 

area  (1.0: 

fori  i  0:  i  ■  total  pts-  1 :  i  -  1 )  •( 

area  area  -  pts  i  0  *  pts  i  -  1  1  : 


area  area  •  pts  total  pts-  I  0  *  ptsO  1 

fori  i  0:  i  total  pts-  1 :  i  -  1 )  { 

area  area  -  pts  i  1  '  pts  i  -  1  0  : 


area  area  -  p’s  total  pis-  1  1  pts  0  0  : 

if!  area  0) 
result  0: 
else 

result  1: 
return  ^result): 


HOI  TINE: 


decompose. r 


Decompose  the  main  polygon  and  as  convex  polygons  are  determined 
these  are  passed  through  calls  to  the  1K1>  system  to  be  filled. 


Object  decompose  poly  goni  total  pts.  pts.  old  ptsl 

int  total  pts: 

(  'oord  pt  s  3  ; 

<  oord  old  pt  s  3  : 


boolean  no  point  found,  searcn  incomplete,  adjust  needed.  OK: 
boolean  skit)  Hag.  skip  search: 
int  pt  1 .  pt'J.  pol\  count,  value: 
point  'working  pt r: 

clear  stack!  ); 

get  ptltotal  pts.  pts): 

get  pt  I  total  pt  ".  pts  i: 

searcn  incomplete  TKl  K: 
skip  search  TKl  !.: 
poly  count  t); 


Once  an  initial  two  points  are  taken  from  the  concave  polygon, 
the  process  moves  through  the  main  chain  to  select  out  points 
for  constructing  the  simple  convex  polygons.  Calls  are  made  to 
various  routines  to  check  for  validity  of  pts  that  are  accepted 
into  a  poly  gon 


w  hilrl  search  incomplete)  { 
adjust  needed  TKl  h: 

OK  TKl  K. 
w  hile( OK  )  { 
get  pt  (total  pts.  pts): 

*  we  first  check  to  see  if  we  have  made  a  complete  loop  of  the 
vertices. 

iff  loop  completedl))  [ 
working  ptr  sirtrt  pt  r : 
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'  find  the  structure  in  the  linked  list  that  matches  the 
vertex  on  top  of  the  stark 

while) working  ptr-  index  !  stark  TOP  .index)  { 
working  ptr  working  ptr-  forward: 

I 

( 

check  to  see  if  vertex  creates  a  line  segment  with  the 
vertex  below  it.  on  the  stark,  that  intersects  any  edges 
of  the  concave  polygon 

iffstack  TOP- 1  .index  !  working  ptr-  back-  index)  ) 
value  intersecting  lines) pts  stark  TOP  .index  0. 

pts  stack  TOP  .index  I  .  pis  stack  TOP-1  .index  0  . 
pts  -,tark  TOP- 1  index  1  .  stack  TOP  .slope, 
stac  k  TOP  .intercept,  total  pts.  pts): 

ill  value)  , 

’  if  vertex  creates  an  intersecting  line  segment  then  we 
remove  it  from  the  stack  and  reset  the  pointer  to  the 
linked  list 

current  [it  r  start  ptr: 

w  hile(  current  ptr-  index  '  stack  TOP- 1  .index)  { 
current  ptr  current  ptr-  forward: 

current  pt  r  current  ptr-  forward: 
release  pt)): 
release  pt  ( ) :  ) 

else  | 

since  the  line  segment  create  is  OK.  we  simply  remove  the 
vertex  from  the  stack  since  it  is  the  same  as  the  one  on 
the  bottom  of  the  stack,  namclv  the  beginning  vertex  ' 

OK  KALSK: 

release  pt  (  ) ; 

if  we  have  not  completed  the  loop  then  we  must  conduct  the 
tests  to  ensure  valid  selection  of  vertices  for  a  convex 


polygon 


else  if(pt  out  of  bounds(totai  pts.  pts))  { 
release  pt(): 

\ 

} 


4  if  at  least  three  points  have  been  selected  then...  * 
i f ( valid  ptsQ)  { 

save  the  beginning  and  ending  vertices  for  future  reference  * 

ptl=  stack  BOTTOM  .index: 
pt2  stack  TOP  .index: 

'  pass  the  vertices  to  IRIS  system  hardware  for  filling  and  display  * 

define  polygon  (old  pts): 
polv  count  -  I: 

4  if  this  is  a  normal  decomposition  and  not  a  "fill  in"  polygon  to 
cover  holes  left  in  the  concave  polygon  the  the  do  the  following 

if(skip  search )  { 

4  set  the  Hags  that  indicate  next  time  through  this  loop  we 
will  be  working  on  a  "fill  in"  polygon  * 

adjust  needed  FALSE: 
skip  search—  FALSE: 

4  adjust  the  pointers  to  the  linked  list  so  next  time  through 
the  loop  we  have  the  start  and  end  vertices  of  the  previous 
polygon  to  begin  working  with 

working  ptr  start  ptr; 
while) working  ptr-  -index  !  pt2)  { 
working  ptr-  working  ptr-  -forward: 

I 

f 


‘  if  the  start  and  end  vertices  of  the  previous  polygon  are 
adjacent  on  the  concave  polygon  then  forget  about  building 
a  "fill  in"  polygon  * 

if) working  ptr-  •  forward-  /index  !  pt  1 )  { 
while(current  ptr-  index  !  pt  1 )  { 
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current  ptr-  current  ptr-  -forward: 

} 

get  pt(total  pts.  pts): 
current  ptr  -  working  ptr: 
get  pt( total  pts.  pts):  } 
else  { 


ROl'TINE: 


clear  stack. c 


This  routine  simply  resets  the  stack  pointer  to  -1.  In  this 
manner  when  the  slack  has  its  first  item  pushed  on.  the  pointer 
will  accurately  reflect  a  value  of  0.  (BOTTOM  OF  STACK) 


clear  stack!  ) 


TOP  -1; 
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A**************** 


ROUTINE:  get  pt.c 


This  retrieves  the  next  point  in  the  main  chain,  as  identified 
by  the  current  ptr.  and  places  the  point  on  the  top  of  the  stack 
In  doing  so  it  increments  the  TOP  variable.  Also  as  the  point 
is  retrieved,  the  slope,  intercept,  and  direction  of  the  line 
segment  it  creates  with  the  vertex  below  it  on  the  stack  is 
calculated. 


get  ptftotal  pts.  pts) 

int  total  pts: 

Coord  pts  3  : 

I 

TOP  1: 

if(TOP  -  total  pts)  { 

*  select  the  vertex  in  the  linked  list  that  is  currently  pointed  to  " 
stack  TOP  .index  current  ptr-  index: 

and  then  advance  the  list  pointer  ' 
current  ptr  current  ptr-  forward: 
iff  TOP  ()){ 

take  the  x.  y.  and  z  coordinates  for  both  the  vertex  on  top  of 
the  stack  and  the  vertex  just  below  it.  to  calculate  the  slope 
and  y-intercept  of  the  line  segment  they  create  * 
determine  slope) pts  stark  TOP-1  .index  0.  pts  stack  TOP-1  .index  1 
pts  stack  TOP  .index  0  .  pts  stack  TOP  .index  1  . 

Astack  TOP  .slope.  Astack  TOP  .intercept. 

Astack  TOP  .direction); 

i 

I 

else  { 

printf("STACK  overflow  occured  while  retrieving  another  point"): 

} 

ret  urn: 
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ROITINE: 


release  pt.c 


This  routine  removes  the  item  on  top  of  the  stark  l>> 
decrementing  the  top  of  slack  pointer  TOP 


ROUTINE: 


slopes. f 


»:**.»*.******.******.**■*.***************  ********  **+**•*■***■**■•***•*:**  I 

This  works  through  the  pts  array.  and  the  stark  array.  By 
taking  value  of  the  TOP  of  stark  and  the  point  beneath  it. 
the  rorresponding  coordinates  are  obtained  from  the  pts  array 
and  a  slope,  intercept  and  dirertion  are  ralrulated.  These 
values  are  then  stored  in  the  stark. 

*.**»*»»**.*  4*  ********  **4»4.*.A*44*-4*  ******  4  4***4  4  *  4  4  *  4  4  ■  •  4*  ******  4 


determine  slope) xl.  y  1.  x‘2.  y‘2.  slope,  interrept.  dirertion) 

Coord  x  1 .  y  1 .  x2.  >  2: 
float  'slope: 

Coord  'intercept: 
int  'direction: 


(  oord  dx.  dy : 

*  calculate  the  delta  x  and  delta  y  for  the  to  points 

dx  x2-  x  1 : 
dy  y  2-  y  1 ; 

*  if  tfie  line  segment  is  vertical  t fieri  use  the  value  99999.0  to 
indicate  such  (avoiding  division  by  zero  problems) 

iff  dx  0)  { 

'slope  99999.0: 

’intercept  x 2 :  } 
else  { 

'slope  dy  dx: 

'intercept  y2  -  (dy  dx)  ’  x2: 

f 


'  based  upon  the  changes  in  delta  x  and  delta  y  we  determine  if  we  are 
traversing  down  a  slope  or  up  a  slope  * 

if(dy  0)  ( 
iff  cl  x  0) 

'direction  1: 
else 

‘direction  0: 


cn 


***#*#******«*.».**»»«»*** 


*  *  * 


HOI  TINh:  compare. c 


This  routine  determines  if  a  point  lies  to  the  left  or  right  of 
a  line  segment.  This  is  done  bv  using  the  x  and  >  coordinates 
of  the  point  with  t hie  slope,  intercept,  and  direction  of  the 
line  segment  to  create  a  parallel  line  segment.  The  y  intercept 
(or  in  the  case  of  vertical  line,  the  x  intercept)  is  used  in 
conjunction  with  the  direction  of  traversal  along  the  line 
segment  to  find  which  side  of  the  line  segment,  the  point  lies 
on. 


compare  boundrv  I  x  1  .v  1  .slope. intercept  .direction  ) 

(  oord  x I .  v  1 ; 

Moat  slope: 

(  oord  intercept 
int  direction: 

i 


float  test  value: 
int  value: 

if  the  line  segment  is  vertical  then  we  use  the  x-intercept  to 
determine  if  a  vertex  lies  to  the  left  or  right  of  another  line 
segment 

if|  slope  99999.0) 

test  value  xl: 

else 

test  value  vl- slope  '  xl: 

based  upon  traversal  up  or  down  a  sloping  line  we  determine  if 
a  vertex  is  to  t hie  left  or  right  of  the  line  segment 

sw  it rh  (direct  ion )  ( 

easel):  ifjtest  value  intercept) 

value  1: 
else 

value  0; 
break 

case  1:  ifltest  value  intercept) 

value  1. 


01 


t  a  a  a  a  ••*■*  +  **<******** 


#*x  +  +a*-*M*.a 


»  *•  »  •  *+*M».irx**  4*  *  *  *  » 


ROUTINE:  pt  OOB.c 

This  routine  takes  th<  vertex  index  from  the  top  of  the  stack 
and  uses  its  x.  v  coordinates  to  determine  if  the  point  lies 
out  of  bounds  with  respect  to  the  convex  polygon  being 
created. 

pt  out  of  bounds! total  pts.  pts) 

int  total  pts: 

( ,'oord  pt  s  3  ; 


( 'oord  x  1 .  \  1 : 
int  pt  r.  \  alue: 
boolean  not  finished: 
point  ‘working  ptr: 

obtain  the  coordinates  for  the  vertex  on  top  of  the  stack  * 

\  1  pts  stack  TOP  .index  0  : 

y  1  pts  stack  TOP  .index  1  : 

ptr  TOP-  I: 

value  0: 

'  use  the  above  coordinates  to  determine  if  the  vertex  lies  to 
the  left  (out  of  bounds)  with  respect  to  any  line  segment 
created  thus  far  in  the  selection  of  vertices  * 

while  (ptr  BOTTOM  lvalue)  { 
value  compare  boundry( x  1  ,y  1  .stack  ptr  .slope.stack  ptr  .intercept. 

stack  ptr  .direction): 
pt  r-  1 . 


If  point  is  not  out  of  bounds  with  any  line  segments  created 
then  we  check  to  see  if  the  segment  produced  by  this  new  point 
addition  causes  the  start  point  for  the  polygon  to  fall  out  of 
bounds 

*  *  •  t  1  «  •  ■  4  4  •  •»4*44»4»*44*««*»*4 


iff  lvalue )  { 

xl  pts  stack  BOTTOM  .index  0  . 
v  I  pts  stack  BOTTOM  .index  1  : 
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valuer  compare  boundry (xl.y  1  .stack  TOP  .slope. stack  TOP  .intercept 
stack  TOP  .direction); 
if(value)  { 
release  pt(): 

current  ptr  current  ptr-  -back: 
working  ptr  -  start  ptr: 

1 


If  the  newly  created  line  segment  does  not  cause  the  initial  pt 
of  the  developing  polygon  to  fall  out  of  bounds,  then  we  check 
to  see  if  the  iirie  seginent  intersects  an>  line  segments  of  the 
concave  polygon. 


iff  lvalue)  { 

working  ptr^  start  ptr; 

w  hile(  working  ptr-  index  !=  stack  TOP. index)  { 
working  ptr  working  ptr-  forward; 

t 

J 

if(stack  TOP-1  .index  I  working  ptr-  back-  -index)  { 

we  use  the  coordinates  for  the  vertex  at  the  top  of  stack 
and  the  vertex  below  it  to  create  a  line  segment  which  we 
then  check  for  intersection  with  any  of  the  edges  of  the 
original  concave  polygon  * 

value  intersecting  linesfpts  stack  TOP  .index  0  . 

pts  stack  TOP  .index  1  .  pts  stack  TOP-1  .index  0  . 
pts  stack  TOP-1  .index  1  .  stack  TOP  .slope, 
stack  TOP  .intercept,  total  pts.  pts): 

\ 

) 


If  the  newly  created  line  segment  does  not  intersect  any  line 
segments  of  the  concave  polygon  then  we  check  to  see  if  the 
line  segment  of  the  original  polygon,  that  the  top  of  stack 
point  heads,  causes  the  point  below  top  of  stack  to  fall  out. 


if(lvalue)  •[ 

if(current  ptr-  back-  back-  index  I  stack  TOP-1  .index)  ) 
xl  pts  stack  TOP- 1  .index  0: 
yl  pts  stack  TOP-1  .index  1  . 

value  compare  boundry ( x  1  ,y  1  .current  ptr-  back-  slope. 


current  ptr- -back-  intercept, 
current  ptr-  back-  direction): 


If  the  point  is  determined  to  be  out  of  bounds,  then  the 
appropriate  pointers  are  adjusted  to  continue  the  decomposition 
process.  Otherwise  the  pt  is  a  valid  addition  to  the  stack 

••*•*•**+**  1**»4**t»**-*-*.»**  *•*•****9*4*4  4*».*»4**»*»*» 

if(value)  { 

working  ptr  start  ptr: 

not  finished  TRIE: 

do  j 

if(working  ptr-  index  stack  TOP. index)  { 

not  finished  FALSE: 
ifjworkirig  ptr-  failed  FALSE) 

working  ptr-  failed  TRI  E:) 
else  j 

working  ptr  working  ptr-  forward: 

I 

f 

}  w hilef not  finished): 


retumf  value): 


ROl'TINE: 


loop  complet.c 


This  routine  checks  to  see  if  tne  vertices  of  the  concave 
polygon  have  been  completely  traversed  vvhiie  developing  a  single 
convex  polygon  This  is  done  when  the  vertex  on  top  of  the 
stack  i>  tne  same  vertex  on  t tie  bottom  of  the  stack.  W  hen  this 
occurs  a  value  of  1  is  returned 


mi  loop  completed!  i 

ir.t  value: 
if:  TOP  -Hi 

t  stark  1  Ol>  index  s;acK  BOTTOM  .index)  • 
value  1; 
else  i 
value  II: 

else  ■ 

print  f(  "hrror  in  stack  index  while  checking  TOP  BOTTOM  values" 
value  -1: 

ret  urn  i  v  alue  I: 
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***«.A**-A*» 


*  *  *  *  *  '0 


*  *  *  * 


ROUTINE:  adjust  .r 

'Phis  routine  searches  t  he  circular  list  of  vert  ices,  checking 
the  various  condition  Hags,  to  locate  a  vertex  or  line  segment 
of  the  originial  polygon  which  has  not  been  utilized.  As  one  is 
located,  t  fit-  list's  currency  pointer  is  moved  to  that  structure 
in  the  list  so  that  further  polygon  decomposition  may  continue 

adjust  search  areal  pis) 

(  oord  pts  .'{  : 


boolean  done: 
point  'pi: 

done  0: 
pi  start  pt r : 

do  •[ 

if(  p  1  -  failed  1  ! 

'  check  to  see  if  the  line  segment  to  either  side  of  the  vertex 
has  been  used  in  a  previously  developed  polygon, 
iflllpl-  forward-  segment  used)  !  ( p  I  -  ■  segment  used))  { 
current  pt  r  pi: 
pi-  •  failed  0: 
done  I:  ) 
else  { 

pi  pi-  •  forward: 

i 

I 

else  { 

pi  pi-  forward: 

I 

i 

\  whilefMone  KK  pi  !  start  ptr): 
ret  urn: 

I 

» 


*X*.*4  +  n4m*.X**  ***■***$*•*.  *»*•***•***»*****  **««!ilt><MlU 

ROl'TINE:  valid  pts.c 

*  •**#****#**#*********#»**#%*»***•***  *».#4*.***4r*«***»**  •  * 

This  routine  checks  to  ensure  that  at  least  3  vertices  have 
been  determined  for  the  convex  polygon,  if  not  then  the 
value  of  0  is  returned. 


boolean  valid  pts( ) 

J 

t 

int  value: 

iff  TOP  2) 
value  1: 
else 

value  0: 
ret  urn  I  value): 
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#*******.*»********♦****  »******»•*♦*****#-♦***»  ».*»«•**•**#**»•*** 
ROl’TINE:  define  polv.r 

*********************************  »  •********♦*»  *  »  »  *  »*«♦**  *  *  »  *  *  • 

This  routine  takes  the  vertices  accumulated  on  the  stack  and 
uses  the  IRIS  features  pmv  and  pdr  to  actually  perform  the 
display  of  the  polygon  created.  The  index  for  each  vertex 
is  used  to  access  the  original  array  of  points  to  get  the  x.  \ 
and  z  coordinates  used. 

*•.*«*.*«••**«•**»««#•  444*****»*****t*-******<<r*»»»»***»**** 


define  polygon (pts) 

Coord  pts  If  : 

< 

ini  count. 

point  'track  ptr: 

track  ptr  start  ptr: 

'  pass  the  starting  vertex  for  t fit-  newly  created  convex 
polygon  to  the  IRIS  hardware  with  a  point  move  command 

pmv(pts  stack  BOTTOM  .index  0  . 
pts  stack  BOTTOM  index  1  . 
pts  stack  BOTTOM  index  2  ): 

'  locate  this  vertex  in  the  linked  list 

whileftrack  ptr-  -index  !-  stack  BOTTOM  index)  { 
track  ptr  track  ptr-  -forward: 

} 


'  set  the  Hag  to  indicate  the  vertex  has  been  used 
in  a  convex  polygon 

iff  track  ptr-  -back-  -index  stark  TOR  .index) 
track  ptr-  -segment  used  TRl’E: 

’  for  each  vertex  of  the  convex  polygon,  pass  its 
coordinates  to  IRIS  hardware  with  the  point  draw 
command 

count  1: 

whilefcount  •  TOR)  { 

whileftrack  ptr-  -index  !  stack  count  .index)  { 
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4**4*** ************ *********«****«»4* ************** 


ROUTINE:  obj  complet.r 

4  **********4*******».*4*******4*4  *  ********  *  *  ****  ***** 

This  routine  checks  each  vertex  and  edge  of  the  original 
concave  polygon  to  see  if  they  have  all  been  utilized  in  the 
process  of  decomposition.  If  so  a  value  of  1  is  returned. 

************  *******  *  *  *  *  **********  *«*4**»*** ***»«**»* 

int  object  complete! total  pts) 
int  total  pts: 


point  '  pi: 

int  count,  value: 

count  0: 
value  1: 
pi  start  ptr: 

search  through  the  linked  list  of  structures,  for  the 
vertices,  and  return  the  index  of  the  first  vertex  that 
creates  an  edge  of  the  concave  polygon  not  yet  utilized  ’ 
do  { 

iff  pi-  ■  segment  used)  j 
pi  pi-  •  forward: 
count  -  1:  } 

else  •) 

value  0: 

I 

I 

)  whilefcount  •  total  pts  value  ~~  1); 
return  (value): 


} 
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ROUTINE: 


intersecting.c 


This  routine  takes  two  vertices  and  the  slope  and  intercept 
of  the  line  segment  they  create  to  determine  if  it  intersects 
any  of  the  edges  of  the  original  polygon.  If  so  then  a  value 
of  1  is  ret  urned. 


intersecting  Iines(x2.  >2.  xl.  y  1.  slope,  intercept,  total  pts.  pts) 

(  oord  \2.  y  2.  x  1 .  y  1 : 
float  slope: 

( ioord  intercept : 
int  total  pts: 

C.oord  pts  3  : 


point  'pi: 

float  x  intersect,  v  intersect,  x  min.  x  max.  \  min.  y  max.  fudge: 
int  value: 

define  a  value  to  offset  round  off  error  ' 

fudge  0.000(H)  1; 

value  0: 
pi  start  ptr: 

for  each  edge  of  the  concave  polygon  we  check  if  either  vertex 
lies  on  the  edge.  If  not  then  we  determine  the  intersection 
of  the  line,  defined  by  the  edge,  with  the  line  defined  by  the 
two  vertices,  f  inally  we  determine  if  this  intersection  lies  on 
both  line  segments.  If  so  then  we  have  found  an  actual  intersection. 


check  if  first  vertex  lies  on  the  polygon’s  edge  being  examined 

i f ( ( ( x  1  !  pts  pi-  -index  0  )  (v  1  !  pts  pi-  -index  1  ))  44 
((xl  !-  pts  pi-  back-  -index  0)  (y  1  !  pts  pi-  -back-  -index  1  ))) 

check  if  second  vertex  lies  on  the  edge 


if(((x2  !-  pts  pi-  -index  0)  ( > 2  !  pts  pi-  index  1))44 

( ( x 2  !  pts  pi-  back-  index  0  )  (>2  ■  pts  pi-  -back-  index  1 


if  the  lines  are  not  parallel  then  ...  ' 
iffslope  !-  pi- aslope)  { 

if  the  line  created  by  the  vertices  is  vertical 
if(s!ope  -  99999.0)  \ 

determine  the  coordinates  of  the  intersection  ‘ 
x  intersect  =  x‘2: 

>  intersect  --  pi-  -slope  '  x  intersect  -  pi-  intercept:  } 

if  the  edge  of  the  polygon  is  vertical  * 

else  i f(  pi-  slope  99999.0)  { 

determine  the  coordinates  of  the  intersection  * 
x  intersect  pts  t>|-  index  0: 

>  intersect  slope  ‘  x  intersect  -  intercept:  ] 

if  neither  line  is  vertical  ’ 
else  { 

determine  the  coordinates  of  the  intersection  ‘ 
x  intersect  (intercept  -  (pi-  interceut))  ((pi-  slope)  -  slope): 
v  intersect  slope  *  x  intersect  -  intercept: 

i 

determine  the  x  range  that  the  intersection  must  occur  iti 
to  lie  on  the  line  segment  created  by  the  two  vertices  * 

max  min(x2.xl.&x  rnax.&rx  min); 

check  if  the  intersection  falls  within  this  range  ‘ 

if((x  intersect  •  x  min-fudge)&&’( x  intersect  ■  x  max -fudge))  { 

since  the  intersection  is  within  this  x  range  we  now  determine 
the  x  range  for  the  edge  of  the  concave  polygon 

max  min(pts  pi-  index  0  .  pts  pi-  back-  index  0  .  k\  max.  k\  min 

check  if  the  intersection  falls  within  this  range  ‘ 

if((x  intersect  x  min-fudgel&’&l  x  intersect  •  x  max -fudge))  { 

now  that  the  intersection  falls  within  the  x  range  of  both 
line  segments,  we  must  calculate  the  \  range  for  each 


segment  and  perform  the  checks 

‘  >  range  for  the  vertices  line  segment 

max  min(y2.y  l.A  y  max. A  s  min): 

if((>  intersect  >  min-fudge)  A  A  ( \  intersect  y  max  fudge))  \ 

>  range  for  t  fie  concave  polygon  edge 

max  minlptspl-  index  1  .pts  pi-  hack-  index  I  ,A>  max. A  \  mm) 
if((>  intersect  >  min-fudge)AA|>  intersect  v  max  -  fudge) ) 
value  I: 

’  the  intersect  ion  exists  for  the  actual  line  segments  and 
therefore  the  vertex  on  top  of  the  stack  must  he  removed  ‘ 


\ 

i 


pi  pi-  forward: 

)while((pl  !  start  ptr)AA(  value  0)): 


ret  urn!  v  alue): 


ROUTINE: 


max  min.c 


Given  two  float  values,  pi  and  p2.  this  routine  determines 
which  is  larger  of  the  two. 

•  *»**»*««***««**4.«*«***«*4****.***»**»***»***»***»**« 

max  rnin(pl.  p2.  max.  min) 

Coord  pi.  p2.  'max.  min : 


iflpl  p'^)  •! 

'max  p2: 
'min  pi:  } 
else  { 

'max  pi. 
'mm  p2: 

\ 


ret  urn: 
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